Magnetic nature of superconductivity in doped cuprates 
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Within the kinetic energy driven superconducting mechanism, the magnetic nature of cuprate 
superconductors is discussed. It is shown that the superconducting state is controlled by both charge 
carrier gap function and quasiparticle coherent weight. This quasiparticle coherent weight grows 
linearly with the hole doping concentration in the underdoped and optimally doped regimes, and then 
decreases with doping in the overdoped regime, which leads to that the maximal superconducting 
transition temperature occurs around the optimal doping, and then decreases in both underdoped 
and overdoped regimes. Within this framework, we calculate the dynamical spin structure factor of 
cuprate superconductors, and reproduce all main features of inelastic neutron scattering experiments, 
including the energy dependence of the incommensurate magnetic scattering at both low and high 
energies and commensurate resonance at intermediate energy. 
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I. INTRODUCTION 

The interplay between the strong electron correlation 
and superconductivity is one of the most important prob- 
lems raised by the discovery of cuprate superconductors^. 
After intensive investigations over more than a decade, 
it has become clear that the strong electron correlation 
in doped cuprates plays a crucial role not only for the 
unusual normal-state behavior but also for the supercon- 
ducting (SC) mechanism^"'^. The parent compound of 
cuprates superconductors is a Mott insulator with the 
antiferromagnetic (AF) long-range order (AFLRO), then 
changing the carrier concentration by ionic substitution 
or increasing the oxygen content turns these compounds 
into the SC-state leaving the AF short-range correlation 
(AFSRC) still intact"*. As a function of the hole doping 
concentration, the SC transition temperature reaches a 
maximum in the optimal doping, and then decreases in 
both underdoped and overdoped regimes^. Moreover, 
this SC transition temperature is dependence of both 
charge carrier gap parameter and quasiparticle coherent 
weight^, which strongly suggests that the quasiparticle 
coherence plays an important role in superconductivity. 

By virtue of systematic studies using the nuclear mag- 
netic resonance, and muon spin rotation techniques, 
particularly the inelastic neutron scattering, the dop- 
ing and energy dependent magnetic excitations in doped 
cuprates in the SC-state have been well established: (a) 
at low energy, the incommensurate (IC) magnetic scat- 
tering peaks are shifted from the AF wave vector [tTjTt] 
to four points [(1 ± (5)7r,7r] and [tt, (1 ± ^)7r] (in units 
of inverse lattice constant) with i5 as the incommensu- 
rability parameter (b) then with increasing energy 
these IC magnetic scattering peaks are converged on 
the commensurate [7r,7r] resonance peak at intermedi- 
ate energy^'*°~*^; and (c) well above this resonance en- 
ergy, the continuum of magnetic excitations peaked at 
IC positions in the diagonal direction are oberved^^"*^. 
It has been emphasized that the geometry of these IC 
magnetic excitations is two-dimensional*^'*'^. Although 
some of these magnetic properties have been observed 



in the normal-state, these IC magnetic scattering and 
commensurate resonance are the main new feature that 
appears into the SC-state. Moreover, AFSRC coexists 
with the SC-state in the whole SC regime^, and the un- 
usual magnetic excitations at high energy have energies 
greater than the SC pairing energy, are present at the SC 
transition temperature, and have spectral weight far ex- 
ceeding that of the resonance*^'*^. These provide a clear 
link between the charge carrier pairing mechanism and 
magnetic excitations in cuprate superconductors. 

Recently, we*^ have discussed the kinetic energy driven 
SC mechanism in doped cuprates based on the charge- 
spin separation (CSS) fermion-spin theory*^, where the 
dressed holons interact occurring directly through the ki- 
netic energy by exchanging dressed spin excitations, lead- 
ing to a net attractive force between dressed holons, then 
the electron Cooper pairs originating from the dressed 
holon pairing state are due to the charge-spin recombi- 
nation, and their condensation reveals the SC ground- 
state. The SC transition temperature is proportional to 
the hole doping concentration in the underdoped regime. 
However, an obvious weakness is that the SC transition 
temperature is too high, and not suppressed in the over- 
doped regime*^. In this paper, we study the magnetic 
nature of the kinetic energy superconductivity in doped 
cuprates along with this line. A short version of this 
work was published earlier*^. One of our main results 
is that the SC transition temperature is suppressed to 
low temperatures by considering the quasiparticle coher- 
ence, and therefore the SC transition temperature is con- 
trolled by both charge carrier gap function and quasipar- 
ticle coherent weight. This quasiparticle coherent weight 
is closely related to the dressed holon self-energy from 
the dressed spin pair bubble, and grows linearly with in- 
creasing doping in the underdoped and optimally doped 
regimes, then decreases with increasing doping in the 
overdoped regime, which leads to that the maximal SC 
transition temperature occurs around the optimal dop- 
ing, and then decreases in both underdoped and over- 
doped regimes. Within this SC mechanism, we give a 
theoretical explanation of inelastic neutron scattering ex- 
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periments on cupratc superconductors^' in terms of 
the collective mode in the dressed holon particle-particle 
channel. 

The paper is organized as follows. The interplay be- 
tween the quasiparticle coherence and superconductivity 
is discussed in Sec. II. In Sec. Ill, we calculate explic- 
itly the dynamical spin structure factor of cuprate su- 
perconductors, and reproduce all main features found in 
experiments in the SC-statc^"^"'', including the energy de- 
pendence of the IC magnetic scattering at both low and 
high energies and commensurate [tt, tt] resonance at in- 
termediate energy. Sec. IV is devoted to a summary and 
discussions. 



II. INTERPLAY BETWEEN THE 
QUASIPARTICLE COHERENCE AND 
SUPERCONDUCTIVITY 

In doped cuprates, the single common feature is the 
presence of the two-dimensional Cu02 plane^, it is be- 
lieved that the relatively high SC transition temperature 
is closely related to doped Cu02 planes. It has been ar- 
gued that the essential physics of the doped Cu02 plane 
is contained in the t-J model on a square lattice^, 

H = -t ^ Cl^Ci+fia + M X/ ^io-C'iCT + J ^ Si • Sj+f), 
if}(7 icr irj 

(1) 

with fj = ±x, ±y, (Cicr) is the electron creation (an- 
nihilation) operator, Sj = CjaCi/2 is spin operator with 
^ = {o'x, (^yi<^z) as Pauli matrices, and /tx is the chemical 
potential. The t-J model (1) is subject to an impor- 
tant local constraint to avoid the double occupancy, i.e., 
So- ^ia^io- ^ 1- III model, the strong electron 

correlation manifests itself by this single occupancy local 
constraint, and therefore the crucial requirement is to 
impose this local constraint. This local constraint can be 
treated properly in analytical calculations within the CSS 
fcrmion-spin theory where the constrained electron op- 
erators are decoupled as, = h\-^ Sr and = 

with the spinful fermion operator hia = e~'**''/ii de- 
scribes the charge degree of frciodom together with some 
effects of the spin configuration rearrangements due to 
the presence of the hole itself (dressed holon), while the 
spin operator Si describes the spin degree of freedom 
(dressed spin), then the electron local constraint for the 
single occupancy is satisfied in analytical calculations^^. 
In this CSS fermion-spin representation, the low-energy 
behavior of the t-J model (1) can be expressed as^^ 

H = -tJ2ih.,S+hl^^S7^^ + h,,S7hl-^St^-) 

i'i) 

- + JeS^^i -^i+r), (2) 

i(T ifj 



with Jeff = (1 - .x)2j, and x = {h\^hi„) = {h]hi) is 
the hole doping concentration. As a consequence, the 
kinetic energy (t) term in the t-J model has been ex- 
pressed as the dressed holon-spin interaction, which re- 
flects that even kinetic energy term in the t-J model has 
strong Coulombic contributions due to the restriction of 
single occupancy of a given site. This dressed holon- 
spin interaction is quite strong, and we^^'^^ have shown 
in terms of Eliashberg's strong coupling thcory^°'^^ that 
in the case without AFLRO, this interaction can induce 
the dressed holon pairing state (then the electron Cooper 
pairing state) by exchanging dressed spin excitations in 
the higher power of the hole doping concentration x. 
The angle resolved photoemission spectroscopy (ARPES) 
measurements^^ have shown that in the real space the 
gap function and pairing force have a range of one lattice 
spacing, this indicates that the order parameter for the 
electron Cooper pair can be expressed as, 

= {hi^hr+fiiS^ Sl_^~ - hiihi+ijiS^S^_^_~) 
= -(5+5-+^)Aft, (3) 
with the dressed holon pairing order parameter, 

Afc = {hi+fjihi-f - hi+fj-fhii), (4) 

which shows that the SC order parameter is closely re- 
lated to the dressed holon pairing amplitude, and is 
proportional to the number of doped holes, and not 
to the number of electrons. Following our previous 
discussions^^'^®, the self-consistent equations that satis- 
fied by the full dressed holon diagonal and off-diagonal 
Green's functions are obtained as, 

g{k) = g^'Hk) + 9^'Hm^i\k)9{k) 

-S('^)(-fc)9t(fc)], (5a) 

+ i:i''\-k)g{k)], (5b) 

respectively, where the four- vector notation k = (k, iw„), 
and the dressed holon self-energies arc obtained as, 

x^Y.D('^\p')D(°\p'+p), (6a) 

P,p' iPm 

X i^i)(0)(p')25W(p'+p), (6b) 

where p = (p, ipm), p' = (p', ip'm)^ ^^i*^ ^^ic dressed holon 
and spin mean-field (MF) Green's functions are evaluated 
asi7-i9. 
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P 



(7a) 
(7b) 



with Bp = \[2x,{np - 1) + X(7p - e)], A = 2ZJeff, 
e = 1 + 2*0/ Jeff, 7p = (l/-^)Ei7e'P'^, is the number 
of the nearest neighbor sites, the dressed spin correlation 

functions x = i^i' S^_^-) and Xz = {SfSf^fj), and the MF 
dressed holon and spin excitation spectra are given by, 

= Ztxjk - IJ; (8a) 

+ ^e(^2 - ^axz - Q!X7p)(e - 7p)], (8b) 

where = + (1 - a)/(4Z), A2 = aC + 

(1 — a)/{2Z), the dressed holon particle-hole param- 
eter (j) = {hl^hi+ficr) , and the dressed spin correla- 
tion functions C = (1/^^) E^,v (^iri^r+,') = 
(l/Z2)E._^-,(5f+.5^^^-,). In order to satisfy the sum 

rule of the dressed spin correlation function (S^S^) = 
1/2 in the case without AFLRO, the important de- 
coupling parameter a has been introduced in the MF 
calculation^^'^'*. In the calculation of the sclf-cncrgies 
(6), the dressed spin part has been limited to the MF 
level^^, i.e., the full dressed spin Green's function in Eq. 
(6) has been replaced by the MF dressed spin Green's 
function (7b), since the normal-state charge transport 

obtained at this level can well describe the experimental 

datai8'25. 

Since the pairing force and dressed holon gap func- 
tion have been incorporated into the self-energy function 
'^2^\f^)i then it is called as the eflFective dressed holon 
gap function. On the other hand, the self-energy func- 
tion T,^^\k) renormalizes the MF dressed holon spec- 
trum, and therefore it describes the quasiparticle co- 
herence. Moreover, 's'^\k) is an even function of 
while T,''^\k) is not. In this case, it is convenient to 
break T,^\k) up into its symmetric and antisymmetric 



parts as, E[^\k) = (k) + ioj^^^l' (k) , where T,\l>{k) 
and 'S^(k) are both even functions of itOn- Now we 
define the charge carrier quasiparticle coherent weight 
Zp^{k) = 1 — S^Q^(fc), then the dressed holon diagonal 
and off-diagonal Green's functions in Eq. (5) can be ex- 
pressed as. 



(h), 



9{k) 



iu;n/ZFik)+^k + E[':^{k) 



'^\k) 



(9a) 



(9b) 



As in the conventional superconductor^", the retarded 
function ReS^g^(fc) may be a constant, independent of 



(k, oj). It just renormalizes the chemical potential, and 
therefore can be neglected. Furthermore, we only study 
the static limit of the effective dressed holon gap function 
and quasiparticle coherent weight, i.e., Yi2\k) = A/,(k), 

and Zp^iyi) = 1 — E^Q^(k). In this case, the dressed 
holon diagonal and off-diagonal Green's functions in Eq. 
(9) can be rewritten explicitly as. 



^F(k) 

iujn - -Bk 
^F(k) 
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2 

1 

w„, -I- E\, 



Zf{\^) 



lU)n 



E^ 



(10a) 



(10b) 



with a = 2F(k)ek, Ahz{k) = Zf(k)A^(k), 
and the dressed holon quasiparticle spectrum £k = 

^f^-h I AhzO^) p, this Zplk) reduces the dressed holon 
quasiparticle bandwidth. Although Zp{k) is still a func- 
tion of k, the wave vector dependence is unimportant, 
since everything happens at the electron Fermi surface 
(EFS). In this case, we will approximate ZpCk) by a 
constant, Zp = Zp{\io), where the special wave vector 
ko is defined below. In the CSS fermion-spin theory, 
the electron diagonal Green's function G{i — j,t — t') = 
{{Cia{t);Cj^{t'))) is a convolution of the dressed spin 
Green's function and dressed holon diagonal Green's 
function, which reflects the charge-spin recombination^, 
and in the present case, it can be calculated in terms of 
Eqs. (7b) and (10a) as^^, 
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Gik) = - 



N Z^ 
p 



duj' duj'' 



, 2n 2tt 

np{uj") + nsiui') 
itOn + to" — w' 



As(p,a;')A(p-k,a;" 



^i:MP-k)£;.{(i 



^p — k 



E. 



p-k 



Li(k,p) 



+ 1 
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£^p-k 



^p— k 



+ 



i2(k,p) 



Er. 



- Wr 



i^p-k 

ii(k,p 



i2(k,p) 



lUJr, 



E, 



p-k 



(11) 



where the MF dressed spin spectral function Ag (k, w) = 
— 2ImD(°)(k, cj), the dressed holon spectral function 
A;,(k,cj) = -2Ini9(k,cj), Li(k,p) = np{Ep_^) + 
nsiuip), i2(k,p) = 1 -np{Ep_k.) +nB{u)p), and nsiu)) 
and np{ui) are the boson and fermion distribution func- 
tions, respectively. Then the electron quasiparticle dis- 
persion is determined by the poles of the electron diag- 
onal Green's function (11). At the half- filling, the t-J 
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model is reduced to the AF Heiscnberg model, where 
there is no charge degree of freedom, and then the dressed 
holon excitation spectrum disappears, while the electron 
quasiparticle dispersion is reduced as the spin excitation 
spectrum^^. This electron diagonal Green's function can 
be used to extract the electron momentum distribution 
(then EFS) as^^, 



2 -]vI^"^(P)/_ 2^MP-'i^,^)nF{uj), (12) 



with ns(p) = J^^dujAs{p,uj)ns{Lj)/2'K is the dressed 
spin momentum distribution. Then this electron momen- 
tum distribution can be evaluated explicitly in terms of 
the MF dressed spin Green's function (7b) and dressed 
holon diagonal Green's function (10a) as, 



2 2N ■ 

p 

X ( 1 - |ezlLtanh[i/3i;p_k] 



(13) 



with ni°^(p) = BpC0th(/3a;p/2)/(2wp). Since the dressed 
spins center around [±7r, ±7r] in the Brillouin zone at 
the MF level^^, then the electron momentum distribu- 
tion (13) can be approximately reduced as, 

nk«l/2-pi")Zf(kA-k) 

X [1 - ,ekA-ktanh(/3£;kA_k/2)/i;k^_k]/2, (14) 



with kA = [n,7T], and = (1/A^) Ep=(±.,±.) (p)- 
It has been shown from ARPES experiments^^ that EFS 
is small pockets around [7r/2,7r/2] at small doping, and 
becomes a large EFS at large doping. Therefore in the 
present case the Fermi wave vector from above electron 
momentum distribution can be estimated qualitatively^^ 
as kp ~ [(1 — x)tt/2, (1 — x)tt/2], and is evolution with 
doping. Then the wave vector ko is obtained as ko = 
kA — kp, and we only need to calculate Zp = ZpO^o) as 
mentioned above. Since the charge-spin recombination 
from the convolution of the dressed spin Green's func- 
tion and dressed holon diagonal Green's function leads 
to form EFS^, then the dressed holon quasiparticle co- 
herence Zp appearing in the electron momentum distri- 
bution also reflects the electron quasiparticle coherence. 
We emphasize that the Fermi wave vector kp estimated 
in the present case only is qualitative correct, while the 
quantitative correct EFS obtained within the t-J model 
is an rather complicated problem, and one may need to 
consider the vertex corrections. This and related theo- 
retical ARPES results are under investigations now. 

Some experiments seem consistent with an s-wave 
pairing^*, while other measurements gave the evidence in 
favor of the d-wave pairing^^'"*". This reflects a fact that 
the d-wave gap function oc — belongs to the same 
representation Fi of the orthorhombic crystal group as 



does s-wavc gap function oc + k^. Within some strong 
correlated models, the earlier numerical simulations^^ 
have shown that the s-wave channel was competitive 
with the d-wave, which indicates that superconductiv- 
ity with both s-wave and d-wave symmetries may arise 
directly from the repulsive interactions. For understand- 
ing of these experimental results, we consider both s- 



wave and d-wave cases, i.e., A^^lCk) = A^27k*^ with 
7^*^ = 7k = {coskx+cosky)/2, for the s-wave pairing, and 
Ai'l(k) = \ with 7^'') = (cosA:, - cosfc^)/2, for 

the d-wave pairing, respectively. In this case, the dressed 
holon effective gap parameter and quasiparticle coherent 
weight in Eq. (6) satisfy the following two equations^'''^^, 

1 - (Zt)^— V -y? -Y^"^ (g) B^Bp 

- ^^^> m 2^ Tk+q7k-p+q7k ^^^^ 



7V3 

k,q,p 

i^i^'^(k,q,p) 



+ ■ 



i^f^(k,q,p) 



(Wp - Wq)2 - E\ (Wp Wq 



El 



(15a) 



Z-^ = X^{Ztf—^Y.^l^^oZF 



q,p 



BqBp 



{Wp-Uq - Sp-q+ko)^ 



f^f (q,p) 



+ 



(a;p-Wq + £;p-q+ko)^ 

i^f (q,p) 

(Wp + Wq-iJp-q+ko)^ 

F^'\ci,p) 

(Wp + Wq + iJp-q+ko)^ 



(15b) 



respectively, where 

a = s,d, i^i^^^ (k, q, p) = (tjp - cjq)[nB(wq) - nB(a)p)][l - 
2np{Ei,)] + Ei^[nB{uJp)nB{-ujq) + nB(t^q)"B(-Wp)], 

■F'P'(k,q,p) = -{ujp + Wq)[nB(t^q) - nB{-ojp)][l - 
2np(Ei^)] + Sk[nB(Wp)nB(wq) + riB(-ci,'p)nB(-Wq)], 

^i^'(q:P) = np{Ep^q+i,^)[nB{uJq) - nij(wp)] - 

(2) 

nB(wp)ni3(-Wq), (q,p) = riF (-Bp-q+ko)['^B(wp) - 

"-B(wq)] - nB{iOq)nB{-uJp), F^ (q,P) = 
"'F(-Ep_q+ko)["s(Wq) -ns(-Wp)] +nB(wp)nB(wq), and 

^i''^(q,p) = np{Ep^q+]^„)[nB{-ujq) - UBiujp)] + 
nB(— a>p)nB(— i-iJq). These two equations arc in control of 
the SC order directly, and must be solved simultaneously 
with other self-consistent equations as shown in Ref.^^, 
then all order parameters, decoupling parameter a, and 
chemical potential /i are determined by the self-consistent 
calculation^^. In this case, we obtain the dressed holon 
pair gap function in terms of the off-diagonal Green's 
function (10b) as^^, 
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1 Al"i(k) 



-tanh[-/3£:k], 



(16) 



then the dressed holon pair order parameter in Eq. (4) 
can be evaluated exphcitly as, 



k 



^"^i^^^tanh[i/3i?k]. (17) 



-^^17,19 ]^ave shown that this dressed holon pairing state 
originating from the kinetic energy term by exchang- 
ing dressed spin excitations can lead to form the elec- 
tron Cooper pairing state, where the SC gap function 
is obtained from the electron off-diagonal Green's func- 
tion P{i - j,t - t') = ((C^(i);Cji(^'))), which is a 
convolution of the dressed spin Green's function and 
dressed holon off-diagonal Green's function^, and in the 
present case can be obtained in terms of the dressed spin 
MF Green's function (7b) and dressed holon off-diagonal 
Green's function (10b) as, 

/t(fc) = l^^^^^ri(P-k) i?p 



p k 



(^p -I- -Ep-k)[n_B(cjp) -I- npi-Ep-u.)] 
(icj„)2 - {ujp + £'p-k)^ 
{ujp - -Ep-k)[?^B(^p) + »F(£^p-k)] 
(iw„)2 - (ujp - i?p-k)^ 



, (18) 



then the SC gap function is obtained from this electron 
off-diagonal Green's function as. 



A(-)(k) = -i^/t(k,^..„) 



B 1 



(19) 



which shows that the symmetry of the electron Cooper 
pair is determined by the symmetry of the dressed holon 
pair^*", and therefore the SC gap function can be writ- 
ten as A'^")(k) = A^''^"f^\ then the SC gap parameter 
in Eq. (3) is evaluated in terms of Eqs. (19) and (17) as 

a'") — —xA^ff'- Since the dressed holon (then electron) 
pairing interaction also is doping dependent, then the 
experimental observed SC gap parameter should be the 
effective SC gap parameter A^^^ ~ —x^^^^- In Fig. 1, we 
plot the effective dressed holon pairing (a) and effective 
SC (b) gap parameters in the s-wave symmetry (solid 
line) and d-wave symmetry (dashed line) as a function 
of the hole doping concentration cc at T = 0.002J and 
t/J = 2.5. For comparison, the experimental result'^^ of 
the upper critical field as a function of the hole doping 
concentration is also shown in Fig. 1(b). In a given dop- 
ing concentration, the upper critical field is defined as the 



critical field that destroys the SC-state at the zero tem- 
perature in the given doping concentration, therefore the 
upper critical field also measures the strength of the bind- 
ing of electrons into Cooper pairs like the effective SC 
gap parameter'^^. In other words, both effective SC gap 
parameter and upper critical field have a similar doping 
dependence'^^ . In this sense, our result is in qualitative 
agreement with the experimental data'^^. In particular, 
the value of A^''' increases with increasing doping in the 
underdoped regime, and reaches a maximum in the opti- 
mal doping Xopt ~ 0.18, then decreases in the overdoped 
regime. 

The present result in Eq. (19) also shows that the 
SC transition temperature Ti""* occurring in the case 
of the SC gap parameter A^") = is identical to the 
dressed holon pair transition temperature occurring in 
the case of the effective dressed holon pairing gap param- 
eter AIj"-* = 0. In correspondence with the SC gap pa- 
rameter, the SC transition temperature Ti"'' as a function 
of the hole doping concentration x in the s-wave (solid 



^ 0.4 




0.00 0.05 0.10 0.15 0.20 
X 

FIG. 1. The effective dressed holon pairing (a) and effective 
superconducting (b) gap parameters in the s-wave symmetry 
(solid line) and d-wave symmetry (dashed line) as a function 
of the hole doping concentration in T = 0.002J and t/J = 2.5. 
Inset: the experimental result of the upper critical field as a 
function of the hole doping concentration taken from Ref. 
[32]. 
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line) and d-wave (dashed line) symmetries for t/J = 2.5 
is plotted in Fig. 2 in comparison with the experimental 
result^ (inset). For the s-wave symmetry, the maximal 

SC transition temperature Tc occurs around a partic- 
ular doping concentration x k, 0.11, and then decreases 
in both lower doped and higher doped regimes. How- 
ever, for the d-wave symmetry, the maximal SC transi- 
tion temperature T^f' occurs around the optimal doping 
concentration Xopt ~ 0.18, and then decreases in both 
underdoped and overdoped regimes. Although the SC 
pairing symmetry is doping dependent, the SC state has 
the d-wave symmetry in a wide range of doping, in qual- 
itative agreement with the experiments'^^"'^'^ . Further- 
more, Tc'*' in the underdoped regime (xi*''* in the lower 



(Tc*'' in the lower doped regime) is set by the hole dop- 
ing concentration^^, this reflects that the density of the 
dressed holons directly determines the superfluid density 
in the underdoped regime for the d-wave case (the lower 
doped regime for the s-wave case). Using an reason- 
ably estimative value of J ~ 800K to 12GGK in doped 
cuprates, the SC transition temperature in the optimal 



doping is T, 



(d) 



0.2 J K, 160K ~ 240K, also in quahtative 



agreement with the experimental data^''^^''^^. 

In the present framework of the kinetic energy driven 
superconductivity, the antisymmetric part of the self- 
energy function S^^''(k) (then Zp) describes the quasi- 
particle coherence, and therefore Zp is closely related to 
the quasiparticle density, while the self-energy function 

{h), 



tration x, and therefore "T^c'^ in the underdoped regime 
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doped regime) is proportional to the hole doping concen- v^) describes the effective dressed holon pairing gap 

function. Since the SC-order is established through an 
emerging quasiparticle^, therefore the SC-order is con- 
trolled by both gap function and quasiparticle coher- 
ence, and is reflected explicitly in the self-consistent equa- 
tions (15a) and (15b). To show this point clearly, we 
plot the quasiparticle coherent weight Zp{Tc) as a func- 
tion of the hole doping concentration x for t/ J = 2.5 
in Fig. 3. As seen from Fig. 3, the doping depen- 
dent behavior of the quasiparticle coherent weight re- 
sembles that of the superfluid density in doped cuprates, 
i.e., Zp grows linearly with the hole doping concentra- 
tion in the underdoped and optimally doped regimes, 
and then decreases with increasing doping in the over- 
doped regime, which leads to that the SC transition 
temperature reaches a maximum in the optimal dop- 
ing, and then decreases in both underdoped and over- 
doped regimes. In comparison with Ref.^^, we there- 
fore find that the quasiparticle coherence plays an im- 
portant role in the kinetic energy driven superconduc- 
tivity of doped cuprates. Since cuprate superconductors 
are highly anisotropic materials, therefore the effective 
SC gap function A('*)(k) = A(^)(cosfc^ -|-cosfcj,)/2 for the 
s-wave symmetry or Af'') (k) = Af"*) {cosK - cosky)/2 for 
the d-wave case is dependent on the momentum. Accord- 
ing to a comparison of the density of states as measured 
by scanning tunnelling microscopy'^^ and ARPES spec- 
tral function^ at [vr, 0] point on identical samples, it has 
been shown that the most contributions of the electronic 
states come from [7r,0] point. In this case, although the 
value of the effective SC gap parameter A'^'*^ (then the 

ratio A'^^^/Tc"'') for the s-wave symmetry is larger than 
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FIG. 2. The superconducting transition temperature as a 
function of the hole doping concentration in the s-wave sym- 
metry (solid line) and d-wave symmetry (dashed line) for 
t/J — 2.5. Inset: the experimental result taken from Ref. 
[5]. 
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FIG. 3. The quasiparticle coherent weight Zf(Tc) as a 
function of the hole doping concentration for t/ J = 2.5. 



these A^''' (then the ratio A^'^^/Tc ) in the d-wave case, 
the system has the SC transition temperature Tc'^'' with 
the d-wave symmetry in a wide range of doping. 



III. DOPING AND ENERGY DEPENDENT 
MAGNETIC EXCITATIONS 

In the CSS fcrmion-spin theory, the AF fluctuation 
is dominated by the scattering of the dressed spins"'^^'^^. 
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Since in the normal-state the dressed spins move in the 
dressed holon backgroimd, therefore the dressed spin self- 
energy (then full dressed spin Green's function) in the 
normal-state has been obtained in terms of the collec- 
tive mode in the dressed holon particle-hole channel^^'^^. 
With the help of this full dressed spin Green's func- 
tion in the normal-state, the IC magnetic scattering 
and integrated spin response of doped cuprates in the 
normal-state have been discussed^'^''^''*, and the results of 
the doping dependence of the incommensurability and 
integrated dynamical spin susceptibility are consistent 
with experimental results in the normal-state'*'*. How- 
ever, in the present SC-state discussed in Sec. II, the 
AF fluctuation has been incorporated into the electron 
off-diagonal Green's function (18) (hence the electron 
Cooper pair) in terms of the dressed spin Green's func- 
tion, therefore there is a coexistence of the electron 
Cooper pair and AFSRC, and then AFSRC can persist 
into superconductivity*^. Moreover, in the SC-state, the 
dressed spins move in the dressed holon pair background. 
In this case, we calculate the dressed spin self-energy 
(then the full dressed spin Green's function) in the SC- 
state in terms of the collective mode in the dressed holon 
particle-particle channel, and then give a theoretical ex- 
planation of the IC magnetic scattering peaks at both 
low and high energies and commensurate resonance peak 
at intermediate energy in the SC-state'^'*'*'^"*^. 

Following our previous discussions for the normal-state 
case**'^*, the full dressed spin Green's functions is ex- 
pressed as. 



D(k,u;) 



£'(o)-i(k,a;)-SW(k,a;) 



(20) 



with the second order spin self-energy S^^^ (k, oj). Within 
the framework of the equation of motion method""^^'^^, 
this self-energy in the SC-state with the d-wave sym- 
metry is obtained from the dressed holon bubble in the 
dressed holon particle-particle channel as, 



1 

'J 



p>p' 



X ■^El)W(p' + A;)iES>t(p)cj(p+y), (21) 



and can be evaluated explicitly in terms of the dressed 
holon off-diagonal Green's function (10b) and dressed 
spin MF Green's function (7b) as, 

E(^)(k,a;) = iZtr^J2i<+P+^ + <-^) 



p.q 

72 Af'') 



i^P(k,p,q) 



y-{Ep-Ep+^ + uj^+:,y 



+ 



+ 



+ 



(Ep+q -Ep+ Wq+k)2 



a;2 _ {Ep + E, 

p(4) 



p-(-q 



^q+k)^ 



f'r(k,p,q) 



a;2 _ {Ep+^ + Ep- Wq+k)^ 



(22) 



where 

Fi^^(k,p,q) = [Ep-Ep+^+uj^+i,){nB{uJci+\,)[nF{Ep)- 
nF{Ep+^)] - nF{Ep+c^)nF{-Ep)], F?^(k,p,q) = 
(iJp+q - Ep+ a;q+k){nB(a;q+k)[nF(-Ep+q) - nF{Ep)] - 
nF{Ep)nF{-Ep+^)}, i?'f)(k,p,q) = {Ep + Ep+^ + 

Wq+k){nB(Wq+k)[ni^(-i?p) - 

nF{Ep+^)] H- nF{-Ep+^)nF{-Ep)}, i^?)(k,p,q) = 

(i?p+i?p+q-Wq+k){nB(t^q+k)[ri_F(--Ep) -n_F(i!;p+q)] - 

nF{Ep+q)nF{Ep)} . With the help of this full dressed 
spin Green's function, we can obtain the dynamical spin 
structure factor in the SC-state with the d-wave symme- 
try as. 



S (k, w) = -2[l + ns(w)]ImZ)(k, uj)=2[l + ns(a;)] 

5gImS](«)(k,u;) 

^ [a;2 - w2 _ BkReE(*)(k, a;)]2 + [BkImE(«)(k, a;)]^ 



(23) 



where ImE^^) (k, u.') and ReE(")(k,cj) are the imaginary 
and real parts of the second order spin self-energy in Eq. 
(18), respectively. 

We are now ready to discuss the doping and energy 
dependent magnetic excitations in the SC-state. In Fig. 
4, we plot the dynamical spin structure factor <5'(k, lo) in 
the {kx, ky) plane at the optimal doping Xopt = 0.18 with 
temperature T = 0.002 J for parameter t/J = 2.5 at en- 
ergy (a) 00 = 0.13J, (b) uj = 0.35J, and (c) uo = 0.65J, 
where the distinct feature is the presence of the IC- 
commensurate-IC transition in the spin fluctuation ge- 
ometry. At low energy, the IC peaks are located at 
[(1 ± (5)/2, 1/2] and [1/2, (1 ± 5)/2\ (hereafter we use the 
units of [27r, 2-k\). However, these IC peaks are energy de- 
pendent, i.e., although these magnetic scattering peaks 
retain the IC pattern at [(l±J)/2, 1/2] and [1/2, (l±5)/2] 
at low energy, the positions of the IC peaks move towards 
[1/2, 1/2] with increasing energy, and then the commen- 
surate [1/2,1/2] resonance peak appears at intermedi- 
ate energy = 0.35J. This anticipated resonance en- 
ergy (jjr = 0.35J f» 35 mev (Ref.^^) is not too far from 
the resonance energy ss 41 mev observed in optimally 
doped YBa2Cu306+j,'^'"'^''~'^2. Furthermore, the IC peaks 
are separated again above the resonance energy, and all 
IC peaks lie on a circle of radius of 5' . The values of 
5' at high energy are different from the corresponding 
values of 6 at low energy. Although some IC satellite 
parallel peaks appear, the main weight of the IC peaks 
is in the diagonal direction. Moreover, the separation at 
high energy gradually increases with increasing energy 
although the peaks have a weaker intensity than those 
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FIG. 4. The dynamical spin structure factor ^(k, u>) in the 
ik^,ky) plane at x^pt = 0.18 with T = 0.002J for t/J = 2.5 
at (a) Lj = 0.13J, (b) to = 0.35J, and (c) uj = 0.65J. 



below the resonance energy. To show this point clearly, 
we plot the evolution of the magnetic scattering peaks 
with energy at Xopt — 0.18 in Fig. 5. For comparison, 
the experimental result^^ of YBa2Cu306+j/ with y = 0.7 
{x Ri 0.12) in the SC-state is shown in the same figure. 
The similar experimental results^^'^^ have also been ob- 
tained for YBa2Cu306+iy with different doping concen- 
trations. Our results show that there is a narrow energy 
range for the resonance peak, and therefore the dispersion 
at high energy is distinctly separated from the low energy 
IC fluctuations. The similar narrow energy range for the 
resonance peak has been observed from experiments^^. 
The present results also show that in contrast to the case 
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FIG. 5. The energy dependence of the position of the mag- 
netic scattering peaks at Xopt = 0.18 and T — 0.002J for 
t/J = 2.5. Inset: the experimental result on YBa2Cu306.85 
in the superconducting-state taken from Ref. [12] 
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FIG. 6. The resonance energy ujr as a function of s — Xopt 
in r = 0.002J for t/J — 2.5. Inset: the experimental result 
taken from Ref. [10]. 
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at low energy, the magnetic excitations at high energy 
disperse almost linearly with energy. Furthermore, the 
resonance energy ujr as a function of doping x — Xopt in 
T = 0.002 J for t/J = 2.5 is plotted in Fig. 6 in compari- 
son with the experimental result^° (inset). It is shown 
that in analogy to the doping dependence of the SC 
transition temperature, the magnetic resonance energy 
Wr increases with increasing doping in the underdoped 
regime, and reaches a maximum in the optimal doping, 
then decreases in the overdoped regime. These mediat- 
ing dressed spin excitations in the SC-state are coupled 
to the conducting dressed holons (then electrons) imder 
the kinetic energy driven SC mechanism^^, and have en- 
ergy greater than the dressed holon pairing energy (then 
Cooper pairing energy). We have also made a series of 
scans for ^(k, w) at different temperatures, and found 
that those unusual magnetic excitations are present near 
the SC transition temperature. Although the simple t-J 




K K K 

4 2 5 

FIG. 7. Function VK(k,tj) in x^pt = 0.18 for t/J = 2.5 with 
T = 0.002J from (a) ki = [(1 - 5)/2, 1/2] via k2 = [1/2, 1/2] 
to k3 = [(1 + <5)/2, 1/2] at a; = 0.13.7 (solid line) and 
LO = 0.35J (dashed line), and (b) k* = [(1 - 6')/2, (1 - d')/2] 
via ka = [1/2, 1/2] to kg = [(l-|-<5')/2, {l+S')/2] at w = 0.35J 
(solid line) and uj = 0.65 J (dashed line). 



model can not be regarded as a comprehensive model for 
the quantitative comparison with cuprate superconduc- 
tors, our these results are in qualitative agreement with 
the major experimental observations of doped cuprates 
in the SC-state^'^o-^^ 

The physical interpretation to the above obtained re- 
sults can be found from the property of the renor- 
malized dressed spin excitation spectrum fi^ = + 
Rci;(")(k,f^k) in Eq. (23). Since both MF dressed spin 
excitation spectrum cj^ and dressed spin self-energy func- 
tion E(")(k, uj) in Eq. (22) are strong doping and energy 
dependent, this leads to that the renormalized dressed 
spin excitation spectrum also is strong doping and en- 
ergy dependent. The dynamical spin structure factor in 
Eq. (23) has a well-defined resonance character, where 
5(k, uj) exhibits peaks when the incoming neutron energy 
CO is equal to the renormalized spin excitation, i.e., 



2 12 



0, 



(24) 



for certain critical wave vectors kc = kc^^ at low en- 
ergy, kc = kc^^ at intermediate energy, and kg = 
kc^-* at high energy, then the weight of these peaks 
is dominated by the inverse of the imaginary part of 
the dressed spin self-energy l/ImS('*)(ki^\a;) at low 
energy, I /lmT,^''^\'ki^\ to) at intermediate energy, and 
l/ImE(*)(kc^'', w) at high energy, respectively. In the 
normal-state^^'^^, the dressed holon energy spectrum has 
one branch ^k, while in the present SC-state, the dressed 
holon quasiparticle spectrum has two branches ±£^k, 
this leads to that the dressed spin self-energy function 
S('*^(k, o)) in Eq. (22) is rather complicated, where there 
are four terms in the right side of Eq. (22). In com- 
parison with the normal-state case^^'^^, the contribution 
for the first and second terms in the right side of the 
dressed spin self-energy (22) comes from the lower band 
—Ek of the dressed holon quasiparticle spectrum like the 
normal-state case, while the contribution for the third 
and fourth terms in the right side of the dressed spin self- 
energy (22) comes from the upper band i?k of the dressed 
holon quasiparticle spectrum. During the above calcula- 
tion, we find that the mode which opens downward and 
gives the IC magnetic scattering at low energy is mainly 
determined by the first and second terms in the right 
side of the dressed spin self-energy (22), while the mode 
which opens upward and gives the IC magnetic scatter- 
ing at high energy is essentially dominated by the third 
and fourth terms in the right side of the dressed spin self- 
energy (22), then two modes meet at the commensurate 
[1/2, 1/2] resonance at intermediate energy. This means 
that within the framework of the kinetic energy driven 
superconductivity, as a result of self-consistent motion of 
the dressed holon pairs and spins, the IC magnetic scat- 
tering at both low and high energies and commensurate 
resonance at intermediate energy are developed. This 
reflects that the low and high energy spin excitations 
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drift away from the AF wave vector, or the zero point of 
WikcUj) is shifted from [1/2,1/2] to = k^^' at low 
energy and kc = kc at high energy. With increasing 
energy from low energy or decreasing energy from high 
energy, the spin excitations move towards to [1/2,1/2], 

i.e., the zero point of VF(kc, to) in kg — kc^-* at low energy 
or kc = kc at high energy turns back to [1/2,1/2], then 
the commensurate [1/2, 1/2] resonance appears at inter- 
mediate energy. To show this point clearly, the function 
W{\si,w) in aJopt = 0.18 for t/J = 2.5 with T = 0.002 J 
from (a) ki = [(1 - 6)/2, 1/2] via ka = [1/2, 1/2] to kg = 
[(1 + (5)/2, 1/2] at w = 0.13 J (solid line) and to = 0.35 J 
(dashed line), and (b) k4 = [(1 - S')/2,{1 - 5')/2] via 
ka = [1/2,1/2] toks = [{l+5')/2, {l+S')/2] atw = 0.35J 
(solid line) and u = 0.65 J (dashed line) is plotted in Fig. 
7, where there is a strong angular dependence with actual 
minima in [(1 - 5)/2, 1/2] and [1/2, (1 - (5)/2], [1/2,1/2], 
and [(1 - S')/2, (1 - S')/2] and [(1 + S')/2, (1 + S')/2] for 
low, intermediate, and high energies, respectively. These 
are exactly positions of the IC peaks at both low and 
high energies and resonance peak at intermediate energy 
determined by the dispersion of very well defined renor- 
malized spin excitations. Since the essential physics is 
dominated by the dressed spin self-energy renormaliza- 
tion due to the dressed holon bubble in the dressed holon 
particle-particle channel, then in this sense the mobile 
dressed holon pairs (then the electron Cooper pairs) are 
the key factor leading to the IC magnetic scattering peaks 
at both low and high energies and commensurate reso- 
nance peak at intermediate energy, i.e., the mechanism of 
the IC magnetic scattering and commensurate resonance 
in the SC-state is most likely related to the motion of 
the dressed holon pairs (then the electron Cooper pairs). 
This is why the position of the IC magnetic scattering 
peaks and commensurate resonance peak in the SC-state 
can be determined in the present study within the t-J 
model under the kinetic energy driven SC mechanism, 
while the dressed spin energy dependence is ascribed 
purely to the self-energy effects which arise from the 
the dressed holon bubble in the dressed holon particle- 
particle channel. Our present result in the SC-state and 
the previous result in the normal-state^*''^* show that the 
IC magnetic scattering at low energy appears in both 
SC- and normal-states, this indicates that the IC mag- 
netic scattering at low energy is not associated with the 
SC-state, which is similar to stripe models^^'^*^^^^, where 
the IC magnetic scattering at low energy is due to the 
formation of magnetic domain lines^^'**^"*^. Since the 
commensurate [1/2,1/2] resonance at intermediate en- 
ergy and IC magnetic scattering at high energy are ab- 
sent from the normal-state^*'^^, then our present result 
also show that the commensurate [1/2, 1/2] resonance at 
intermediate energy and IC magnetic scattering at high 
energy are closely related to the SC-state, which is differ- 
ent from the stripe theory^'''^^'^^, where the linear spin 
wave models based the stripe ground state predict that 
the spin excitations at high energy are nearly symmet- 



ric around [1/2, 1/2] position and disperse almost linearly 
with energy^"''^^'^^, then the commensurate [1/2, 1/2] res- 
onance may represent a characteristic energy defined by 
the size of a stripe domain. Although there are some sub- 
tle differences between our present approach and stripe 
theory, both theories can give the qualitative interpreta- 
tion for all main features of the unusual spin response of 
cuprate superconductors^' 

IV. SUMMARY AND DISCUSSIONS 

In summary, within the framework of the kinetic en- 
ergy driven superconductivity^'', we have discussed the 
magnetic nature of cuprate superconductors. It is shown 
that the SC-state is controlled by both SC gap func- 
tion and quasiparticle coherent weight. This quasipar- 
ticle coherent weight is closely related to the dressed 
holon self-energy from the dressed spin pair bubble, and 
grows linearly with the hole doping concentration in the 
underdoped and optimally doped regimes, and then de- 
creases with doping in the overdoped regime, which leads 
to that the maximal SC transition temperature Tc'^^ oc- 
curs around the optimal doping Xopt ~ 0.18, and then 
decreases in both underdoped and overdoped regimes, in 
qualitative agreement with the experiments"''. Although 
the symmetry of the SC-state is doping dependent, the 
SC-state has the d-wave symmetry in a wide range of 
doping. Within this d-wave SC-state, we have calculated 
the dynamical spin structure factor of cuprate supercon- 
ductors in terms of the collective mode in the dressed 
holon particle-particle channel, and reproduce all main 
features of inelastic neutron scattering experiments in the 
SC-state' including the energy dependence of the 
IC magnetic scattering peaks at both low and high ener- 
gies and commensurate resonance peak at intermediate 
energy. In particular, we have shown that the unusual IC 
magnetic excitations at high energy have energies greater 
than the dressed holon pairing energy (then SC Cooper 
pairing energy), and are present at the SC transition tem- 
perature. 

The t-J model is characterized by a competition be- 
tween the kinetic energy (t) and magnetic energy (J). 
The magnetic energy J favors the magnetic order for 
spins, while the kinetic energy f. favors delocalization of 
holes and tends to destroy the magnetic order. Therefore 
the introduction of the additional second neighbor hop- 
ping t' in the t-J model may be equivalent to increase 
the kinetic energy, and this t' term does not change spin 
configuration because of the same sublattice hopping. In 
this case, we**^ have discussed the effect of the additional 
second neighbor hopping t' on superconductivity within 
the t-t'-J model, and found that the d-wave SC pairing 
correlation is enhanced, while the s-wave SC pairing cor- 
relation is heavily suppressed. 

Superconductivity in cuprates emerges when charge 
carriers, holes or electrons, are doped into Mott 
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insulators**^^^. Both holc-dopcd and electron-doped 
cupratc superconductors have the layered structure of 
the square lattice of the Cu02 plane separated by insu- 
lating layers'''^''. In particular, the symmetry of the SC 
order parameter is common in both case^°''*^, manifest- 
ing that two systems have similar underlying SC mech- 
anism. On the other hand, the strong electron correla- 
tion is common for both hole-doped and electron-doped 
cupratcs, then it is possible that superconductivity in 
electron-doped cuprates is also driven by the kinetic en- 
ergy as in hole-doped case. Within the t-t'-J model, 
we"*^ have discussed this issue, and found that supercon- 
ductivity appears around the optimal doping in electron- 
doped cuprates, and the maximum achievable SC transi- 
tion temperature is lower than hole-doped cuprates due 
to the electron-hole asymmetry. 



ACKNOWLEDGMENTS 

The author would like to thank Dr. Ying Liang, Dr. 
Bin Liu, Dr. Jihong Qin, Professor Y.J. Wang, and Pro- 
fessor H.H. Wen for the helpful discussions. This work 
was supported by the National Natural Science Founda- 
tion of China under Grant Nos. 10125415 and 90403005. 



1 P.W. Anderson, Science 235, 1196 (1987). 

^P.W. Anderson, Phys. Rev. Lett. 67, 2092 (1991); Sci- 
ence 288, 480 (2000); Physica C 341-348, 9 (2000); cond- 
mat/0108522. 

^ R.B. Laughlin, Phys. Rev. Lett. 79, 1726 (1997); J. Low. 

Tern. Phys. 99, 443 (1995). 
■* See, e.g., M.A. Kastner, R.J. Birgeneau, G. Shiran, and Y. 

Endoh, Rev. Mod. Phys. 70, 897 (1998). 
^ See, e.g., J.L. Tallon, J.W. Loram, J.R. Cooper, C. 

Panagopoulos, and C. Bernhard, Phys. Rev. B 68, 180501 

(2003). 

® H. Ding, J.R. Engelbrecht, Z. Wang, J.C. Campuzano, S.C. 

Wang, H.B. Yang, R. Rogan, T. Takahashi, K. Kadowaki, 

and D.G. Hinks, Phys. Rev. Lett. 87, 227001 (2001); R.H. 

He, D.L. Feng, H. Eisaki, J. -I. Shinioyama, K. Kishio, and 

G.D. Gu, Phys. Rev. B 69, 220502 (2004). 

P. Dai, H.A. Mook, R.D. Hunt, and F. Dogan, Phys. Rev. 

B63, 54525 (2001); H. He, P. Bourges, Y. Sidis, C. Ul- 

rich, L.P. Rcguault, S. Pailhos, N.S. Borzigiarova, N.N. 

Kolesnikov, and B. Kcimcr, Science 295, 1045 (2002); N.B. 

Christensen, D.F. McMorrow, H.M. R0nnow, B. Lake, S.M. 

Hayden, G. Acppli, T.G. Perring, M. Mangkorntong, N. 

Nohara, and H. Tagaki, Phys. Rev. Lett. 93, 147002 (2004). 
® K. Yamada, C.H. Lee, K. Kurahashi, J. Wada, S. Waki- 

moto, S. Ueki, H. Kimura, Y. Endoh, S. Hosoya, and G. 

Shirane, Phys. Rev. B 57, 6165 (1998). 
® S. Wakimoto, H. Zhang, K. Yamada, I. Swainson, H. Kim, 

and R.J. Birgeneau, Phys. Rev. Lett. 92, 217004 (2004); 



M. Fujita, K. Yamada, H. Hiraka, P.M. Gchring, S.H. Loo, 
S. Wakimoto, and G. Shirane, Pliys. Rev. B 65, 064505 
(2002); S. Wakimoto, G. Shirane, Y. Endoh, K, Hirota, S. 
Ueki, Y.S. Lee, P.M. Gehring, and S.H. Lee, Phys. Rev. B 
60, R769 (1999). 

^° P. Bourges, B. Keimer, S. Pailhos, L.P. Rcgnault, Y. Sidis, 
and C. Ulrich, Physica C 424, 45 (2005); H. He, Y. Sidis, 
P. Bourges, G.D. Gu, A. Ivanov, N. Koshizuka, B. Liang, 
C.T. Lin, L.P. Regnault, E. Schoenherr, and B. Keimer, 
Phys. Rev. Lett. 86, 1610 (2001). 

" P. Bourges, Y. Sidis, H.F. Fong, L.P. Regnauh, J. Bossy, 
A. Ivanov, and B. Keimer, Science 288, 1234 (2000). 
M. Aral, T. Nishijima, Y. Endoh, T. Egami, S. Tajima, 
K. Tomimoto, Y. Shiohara, M. Takahashi, A. Garret, and 
S.M. Bennington, Phys. Rev. Lett. 83, 608 (1999). 

" S.M. Haydon, H.A. Mook, P. Dai, T.G. Perring, and F. 
Dogan, Science 429, 531 (2004). 

" C. Stock, W.J. Buyers, R.A. Cowley, PS. Clegg, R. Coldea, 
CD. Frost, R. Liang, D. Poets, D. Bonn, W.N. Hardy, and 
R.J. Birgeneau, Piiys. Rev. B71, 24522 (2005). 

" J.M. Tranquada, H. Woo, T.G. Perring, H. Goka, G.D. Gu, 
G. Xu, M. Fujita, and K. Yamada, Nature 429, 534 (2004); 
J.M. Tranquada, cond-mat/0512115. 

V. Hinkov, S. Pailhcs, P. Bourges, Y. Sidis, A. Ivanov, 
A. Kulakov, C.T. Lin, D.P. Chen, C. Bernhard, and B. 
Keimer, Science 430, 650 (2004). 

" Shiping Feng, Phys. Rev. B 68, 184501 (2003). 

Shiping Feng, Jihong Qin, and Tianxing Ma, J. Phys. Con- 
dens. Matter 16, 343 (2004); Shiping Feng, Tianxing Ma, 
and Jihong Qin, Mod. Phys. Lett. B17, 361 (2003). 
Tianxing Ma, Huaiming Guo, and Shiping Feng, Mod. 
Phys. Lett. B18, 895 (2004). 

20G.M. Eliashberg, Sov. Phys. JETP 11, 696 (1960); D.J. 
Scalapino, J.R. Schrieffer, and J.W. Wilkins, Phys. Rev. 
148, 263 (1966). 

See, e.g., G.D. Maham, Many Particle Physics, (Plenum 
Press, New York, 1981), Chapter 9. 

22Z.X. Shen, D.S. Dessau, B.C. Wells, D.M. King, W.E. 
Spicer, A.J. Arko, D. MarshaU, L.W. Lombardo, A. Ka- 
pitulnik, P. Dickinson, S. Doniach, J. DiCarlo, T. Loeser, 
and C.H. Park, Phys. Rev. Lett. 70, 1553 (1993); H. Ding, 
M.R. Norman, J.C. Campuzano, M. Randcria, A.F. Bell- 
man, T. Yokoya, T. Takaiiashi, T. Mocliiku, and K. Kad- 
owaki, Phys. Rev. B54, R9678 (1996). 

23 Shiping Feng and Yun Song, Phys. Rev. B 55, 642 (1997). 

2* J. Kondo and K. Yamaji, Prog. Theor. Phys. 47, 807 
(1972). 

2^ Shiping Feng and Zhongbing Huang, Phys. Lett. A 232, 

293 (1997); Fong Yuan, Jihong Qin, Shiping Feng, and 

W.Y. Chen, Phys. Rev. B 67, 134505 (2003). 
2^ Huaiming Guo and Shiping Feng, cond-mat/0509508. 

Z.X. Shen and D.S. Dessau, Phys. Rep. 253, 2 (1985). 
2* P. Chaudhari and S.Y. Lin, Phys. Rev. Lett. 72, 1084 

(1994); D.H. Wu, J. Mao, S.N. Mao, J.L. Peng, X.X. 

Xi, T. Venkatesan, R.L. Greene, and S.M. Anlago, Phys. 

Rev. Lett. 70, 85 (1993); S.M. Anlage, B.W. Langley, G. 

Deutscher, J. Halbritter, and M.R. Beasley, Phys. Rev. B 

44, 9764 (1991). 
2^ J.A. Martindale, S.E. Barrett, K.E. OHara, CP. SUchter, 

W.C Lee, and D.M. Ginsberg, Phys. Rev. B 47, 9155 



11 



(1993); W.N. Hardy, D.A. Bonn, D.C. Morgan, R. Liang, 
and K. Zhang, Phys. Rev. Lett. 70, 3999 (1994); D.A. Woll- 
man, D.J. Van Harlingen, W.C. Lee, D.M. Ginsberg, and 
A.J. Leggett, Phys. Rev. Lett. 71, 2134 (1993). 
See, e.g., C.C. Tsuei and J.R. Kirtley, Rev. Mod. Phys. 72, 
969 (2000). 

E. Dagotto and J. Riera, Phys. Rev. B 46, 12084 (1992); 
R.T. Scalettar, Physica C 162-164, 313 (1989); R.T. 
Scalettar, S.R. White, D.J. Scalapino, and R. Sugar, Phys. 
Rev. B 44, 770 (1991); M. Capone, M. Fabrizio, C. Castel- 
lani, and T. Tosatti, Science 296, 2364 (2002). 
H.H. Wen, H.P. Yang, S.L. Li, X.H. Zeng, A.A. Soukias- 
sian, W.D. Si, and X.X. Xi, Europhys. Lett. 64, 790 (2003). 
N.-C. Yeh, C.T. Chen, G. Hammerl, J. Mannhart, A. 
Schmehl, C.W. Schneider, R.R. Schulz, S. Tajima, K. 
Yoshida, D. Garrigus, and M. Strasik, Phys. Rev. Lett. 
87, 087003 (2001); G. Deutscher, Nature 397, 410 (1999). 
A. Biswas, P. Fournier, M.M. Qazilbash, V.N. Smolyani- 
nova, H. Balci, and R.L. Greene, Phys. Rev. Lett. 88, 
207004 (2002). 

C.C. Tsuei, J.R. Kirtley, G. Hammerl, J. Mannhart, H. 
Raffy, and Z.Z. Li, Phys. Rev. Lett. 93, 187004 (2004). 
Y.J. Uemura, G.M. Luke, B.J. Sternlieb, J.H. Brewer, J.F. 
Carolan, W.N. Hardy, R. Kadono, J.R. Kcmpton, R.F. 
Kiefl, S.R. Kreitzman, P. Mulhern, T.M. Riscrrian, D.L. 
Williams, B.X. Yang, S. Uchida, H. Takagi, J. Gopalakrish- 
nan, A.W. Sleight, M.A. Subramanian, C.L. Chien, M.Z. 
Cieplak, G. Xiao, V.Y. Lee, B.W. Statt, C.E. Stronach, 
W.J. Kosslcr, and X.H. Yu, Phys. Rev. Lett. 62, 2317 
(1989); Y.J. Uemura, L.P. Le, G.M. Luke, B.J. Sternlieb, 
W.D. Wu, J.H. Brewer, T.M. Riscman, C.L. Seaman, M.B. 
Maple, M. Ishikawa, D.G. Hinks, J.D. Jorgensen, G. Saito, 
and H. Yamochi, Phys. Rev. Lett. 66, 2665 (1991). 
Y. DeWilde, N. Miyakawa, P. Guptasarma, M. lavarone, 
L. Ozyuzer, J.F. Zasadzinski, P. Romano, D.G. Hinks, C. 
Kendziora, G.W. Crabtree, and K.E. Gray, Phys. Rev. 
Lett. 80, 153 (1998). 

Shiping Feng and Zhongbing Huang, Phys. Rev. B 57, 
10328 (1998); Feng Yuan et al., Phys. Rev. B 64, 224505 
(2001). 

S. Shamoto, M. Sato, J.M. Tranquada, B.J. Sternlib, and 

G. Shirane, Phys. Rev. B 48, 13817 (1993). 

J. Zaanen and O. Gunnarsson, Phys. Rev. B 40, 7391 

(1989); D. Poilblanc and T.M. Rice, Phys. Rev. B 39, 9749 

(1989). 

E. W. Carlson, D.X. Yao, and D.K. Campbell, Phys. Rev. 
B 70, 064505 (2004); CD. Batista, G. Ortiz, and A.V. 
Balatsky, Phys. Rev. B 64, 172508 (2001). 

F. Kriiger and S. Scheidl, Phys. Rev. B 67, 134512 (2003). 
Shiping Feng and Tianxing Ma, Phys. Lett. A, to be pub- 
lished, cond-mat/0506114. 

Y. Tokura, H. Takagi, and S. Uchida, Nature 337, 345 
(1989); L. AlfT, Y. Krockenberger, B. Welter, M. Scho- 
necke, R. Gross, D. Manske, and M. Naito, Nature 422, 
698 (2003). 

C.C. Tsuei and J.R. Kirtley, Phys. Rev. Lett. 85, 182 

(2000). 

Tianxing Ma and Shiping Feng, Phys. Lett. A339, 131 
(2005). 



12 



